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ABSTRACT 

Due to orbital decay by gravitational-wave radiation, some close-binary helium white 
dwarfs are expected to merge within a Hubble time. The immediate merger prod- 
ucts are believed to be helium-rich sdO stars, essentially helium main-sequence stars. 
We present new evolution calculations for these post-merger stars beyond the core 
helium-burning phase. The most massive He-sdO's develop a strong helium-burning 
shell and evolve to become helium-rich giants. We include nucleosynthesis calcu- 
lations following the merger of O.4M0 helium white-dwarf pairs with metallicities 
Z = 0.0001,0.004,0.008 and 0.02. The surface chemistries of the resulting giants are 
in partial agreement with the observed abundances of R Coronae Borealis and extreme 
helium stars. Such stars might represent a third, albeit rare, evolution channel for the 
latter, in addition to the CO+He white dwarf merger and the very-late thermal pulse 
channels proposed previously. We confirm a recent suggestion that lithium seen in 
R CrB stars could form naturally during the hot phase of a merger in the presence of 
3 He from the donor white dwarf. 

Key words: stars: peculiar (helium), stars: evolution, stars: white dwarfs, stars: 
abundances, binaries: close 



1 INTRODUCTION 

The existence of close binary white dwarfs has been 
predicted theoretically ( Han 



ture, luminosity and surface composition of compact helium- 
rich subdwarf O stars (Zhang & JefFery 20121. In particu- 



1998 Nelemans et al. 



2000| |2001 1 and demonstrated observati onally jRebassa- 



I Mansergas et al.|[20TT| |Brown et al.pOTTj [Kllic et al.||2010[ 
|2011[ ). It has been demonstrated that the orbits of close bi- 
naries decay through the emission of gravitational wave radi- 
ation ( |Evans et al.|1987| [Cropper et al.|1998[ ). Models for the 
merger of nearly-equal mass white dwarfs following spiral- 
in demonstrate that these occur on a dynamical timescale 
(minutes) ( Benz et al.||1990 Guerrero et al. 2004 Pakmor 
et al.||2011[ ). The less massive WD is disrupted and part of 
the debris forms a prompt hot corona, the remainder forms 
a disk which subsequently accretes onto the surviving white 
dwarf ( |Yoon et al.|2007| [Loren-Aguilar et al.|2009[ ). Models 
for the evolution of the product of a double helium-white- 
dwarf merger demonstrate the off-centre ignition of helium 
burning, followed by expansion and evolution onto the he- 
lium main sequence (Saio & JefFery 20001. Models which 
include the composite nature of the debris (corona + disk) 
successfully account for the distribution in effective tempera- 
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lar, the most-massive mergers become the hottest sdO stars 
with carbon-rich surfaces, whilst the least-massive mergers 
are cooler with predominantly nitrogen-rich surfaces. 

Following inward migration of the helium-burning shell, 
the merged star is essentially a helium main-sequence star 
having a mass between 0.4 and 1.0 M@, although its surface 
layers will carry a record of the merger event. The evolu- 
tion of low-mass helium stars has been studied for over 40 
years JPaczyhski 1971 Dinger 1972 Trimble & Paczynski 
1^73}|Weissfl987] . The lowest -mass stars will evolve directly 



to become CO/He white dwarfs following core helium ex- 
haustion. For M ^ I.OMq, shell helium burning will ignite 
around the carbon/oxygen core, and the star will expand 
to become a giant. Early authors proposed this as a possi- 
ble origin for the hydrogen-deficient R Coronae Borealis (R 
CrB) variables ( |Paczyhs^ki|197l] | WeissfT987 1 . 

R CrB stars and the hotter extreme helium stars are 
low-mass supergiants of spectral types F, A and B. Their 
surfaces are extremely deficient in hydrogen, and enriched 
in carbon, oxygen, neon and other nuclear waste. Two prin- 
ciple evolution channels have been established. A small frac- 
tion may be produced following a late thermal pulse in a 
post-asympotic giant-branch star on the white dwarf cool- 
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ing sequence (Iben & Renzini 1984 Clayton et al. 2011 1. The 



majority are more likely to have been produced following the 
merger of a carbon-oxygen white dwarf with a helium white 



dwarf (|Webbink|p84| |Saio fc Jeffeiy| |2002| |Jeffery et al 
2011||Longland et al.|20TT i 



In this letter we demonstrate that high-mass double 
helium-white-dwarf mergers could contribute a third evo- 
lution channel and account for some, at least, of the lowest 
luminosity R CrB and extreme helium stars. 



2 MODELS 
2.1 Evolution 

Numerical simulations of stellar evolution are carried out us- 
ing the stellar-evolution code "Modules for Experiments in 
Stellar Astrophysics" (MESA) ( |Paxton et al.|2011[ ). Helium- 
white-dwarf models were generated by evolving a 2 Mq main 
sequence star and removing the envelope when this reached 



the required mass (see Zhang & Jeffery (20121 for details). 
The post-merger evolution was modelled by considering fast 
accretion at 10 4 Mgyr -1 (representing the formation of a hot 
corona), followed by slow accretion at 10 -5 Mgyr -1 (repre- 
senting accretion from a Keplerian disk). |Zhang fc Jeffery] 
(20121 considered three sets of assumptions concerning this 
accretion: i) all fast, ii) all slow, and iii) composite. In the 
composite model, most of the mass of the donor white dwarf 
is ingested into the hot corona, with 0.1 Mq remaining in the 
Keplerian disk for slow accretion. 

The evolution from off-center helium ignition to the end 
of core helium burning has been discussed and compared 
favourably with observations of helium-rich subdwarf B and 
O stars ( |Zhang fe Jeffery||2l)l"2| ) . Following core helium ex- 
haustion, the post-merger models with M — 0.5,0.6 and 
0.7 Mq developed a brief phase of thick helium shell burn- 
ing, before helium burning ends and the stars contract to be- 
come CO/He white dwarfs. The models with M = 0.8 M 
(Mc OI .c = 0.64 Mq) developed a thin helium-burning shell 
and evolved to high luminosity and low effective tempera- 
ture. This is a regime occupied by another class of hydrogen- 
deficient star, the RCrB stars. 

We have made additional calculations for these 0.4+0.4 
M mergers with metalicities Z = 0.0001,0.004,0.008 and 
0.02 to investigate whether the He+He WD merger channel 
could produce some of the observed EHe and RCrB stars. 



2.2 Nucleosynthesis 

In the SPH simulations of white dwarf mergers, matter in 
the hot coronae may briefly reach temperatures of 6 x 10 8 K 
or more ftoon et al~1|2007| |Loren-Aguilar et al.||2009| |, so 
that some nucleosynthesis of a-rich material will occur. Sim- 
ilarly, our one-dimensional quasi-equilibrium calculations in- 
dicate coronal temperatures up to 4 x 10 8 K for a 0.4+0.4 
M©merger | |Zhang fc Jeffery[2 012) . At these temperatures, 
the 3a and other alpha capture reactions are ignited at the 
surface of the accretor almost immediately. Thus, carbon is 
produced by helium burning and nitrogen is destroyed by a 



capture. As Warner ( 1967 I and Clayton et al 

^N(a, 7 ) 



cated, the destruction of 4 N by the 
becomes more efficient than it is at low temperature 
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Figure 1. Evolutionary tracks at different metallicities on a sur- 
face gravity - effective temperature diagram for a 0.8 Mq post- 
merger model for a He+He white-dwarf binary. The dashed line 
shows the evolution track during accretion, solid lines show the 
post-merger evolution. The colour from grey to dark (a colour 
version is available in the online journal, from pink to red) shows 
the different metallicities, i.e. Z = 0.0001, 0.004, 0.008, .02. The 
symbols shows the R CrB stars from |Jeffery et al 



the circle symbols shows the EHe stars from | Jeffery et al 



(2011), 



120111 



the dots show the He-sdO stars from Hirsch (2009). The inset 
top left shows an enlargement around the phase of core helium 
burning. The inset lower right shows an enlargement of the R CrB 
star region. Three stages are marked: A for inward shell-burning; 
B for outward shell-burning; C for giant branch to white dwarf. 



the temperature continues to increase, ls O produces 22 Ne 
through a capturing. The question is whether the composi- 
tion of these surface layers also matches that of the RCrB 
stars. 



3 RESULTS 

3.1 Evolution tracks 

Taking the Z — 0.02 model as an example, during the fast 
accretion stage a hot corona with a radius of 0.06 Rq is 
formed within 16 minutes, whilst the luminosity (logL/Lg) 
increases from -2 to -1. As the fast accretion starts, the coro- 
nal temperature reaches 10 8 K very quickly. At the end of the 
fast phase, the helium-burning shell reaches a peak tempera- 
ture of about 4 x 10 8 K. During this process, an almost com- 
pletely convective envelope is produced due to the corona 
being heated by the energy produced by helium buring. 

After the fast merger phase, the remaining mass of the 
secondary forms a Keplerian disc surrounding the accretor 
and transfers mass to the central object at a rate comparable 
to the Eddington accretion rate. During this slow-accretion 
phase, the helium-burning shell still heats the corona, forc- 
ing it to continue expanding. Thus, the star expands to 
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~ 32 Rq within 4 x 10 3 yrs and then contracts away from 
the giant branch (Fig. [T]). 

During subsequent evolution, there follow about 20 he- 
lium flashes in 6 x 10 s yrs, each subsequent flash decreasing 
in intensity, until the helium-burning shell (flame) reaches 
the centre of the star. This corresponds to the end of the 
loops on the g — Teff diagram. After this, a standard core- 
helium burning phase is established, followed by a normal 
helium-shell burning phase, leading to higher luminosity, 
and finally cooling to the white dwarf sequence. Just be- 
fore reaching the white dwarf sequence, there is a small loop 
caused by a weak final helium shell flash (Fig. [TJ). 

In order to compare with observations, data for 13 
RCrB stars and 21 EHe stars are taken from|Jeffery et al.] 



( |2011h and for 15 He-sdO stars from |Hirsch| |2009|. As 
Fig. [T] shows, most of the EHe and RCrB star lie around 
the low-gravity evolution corresponding to outward shell- 
helium burning after core-helium burning. Four or five of 
the EHe stars lie closer to the high-gravity inward shell- 
burning phase. After this phase they evolve to become He- 



sdO stars as Saio fe Jeffery ( 2000 1 indicated. Two RCrB 



stars (RT Nor and RS Tel) have high enough gravities to be 
immediate post-merger objects. They may still have a disc 
surrounding them or be surrounded by dust. Thus, there is a 
channel from the the He+He white dwarf merger which can 
produce RCrB stars, then possibly EHe stars, and then He- 
sdO stars. Subsequently, massive He-sdO stars might evolve 
to become EHe and RCrB stars once again, before finally 
cooling to become white dwarfs. We note that this channel 
will not generate metal-poor RCrB stars, due to lower opac- 
ity in the stellar envelope, but can generate low-metallicity 
EHe stars (e.g. BD+10 2179: Kupfer et al. 2012). 

3.2 Surface composition 

At the beginning of the fast-accretion phase, the corona is 
fully convective, and helium burning makes the corona rich 
in 12 C, 18 and 22 Ne. At the end of fast accretion and 
beginning of slow accretion, flash-driven convection mixes 
12 C throughout the 0.3Mq envelope. Although this mate- 
rial is subsequently buried by C-poor material from the disk, 
the deep opacity-driven convection which follows each shell- 
pulse dredges carbon-rich material to the new surface, so 
that the final product, after slow-accretion terminates, is 
rich in 12 C and 22 Ne. 

The most recent measurements of the surface abun- 
dances of EHe and R CrB stars have been collated by |Jef-| 



fery et al. ( 2011 1. Here we summarise the main observational 
features and compare these with our results; see Fig. [2] for 
details. 

Carbon is enriched in all EHe and RCrB stars. Exclud- 
ing MVSgr, the EHes shows a mean carbon abundance 
loge=9.3, and a range from 8.9-9.7. The carbon abundance is 
more difficult to measure in RCrBs. Furthermore, in RCrB 
stars cool enough to show CO, the 12 C/ 13 C ratios are very 
large (500), indicating a 3a or helium-burning origin for the 
carbon excess. In our models, the 12 C was produced by he- 
lium burning through 3a reaction during the hot accretion 
phase, and then brought up to the surface during the slow- 
accretion phase by flash-driven convection. The 13 C abun- 
dance remained almost unchanged during the whole simu- 
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Figure 2. The observed surface abundances of RCrB and EHe 
stars jjeffery et al.|2011} compared with our nucleosynthesis com- 
putation. The axes [X] and [Fe] give logarithmic abundances rel- 
ative to solar for individual elements and for iron, respectively. 
Solid lines are abundances given by our simulation, which include 
four metallicities z = 0.0001,0.004,0.008 and 0.02. The asterisks 
correspond to the solar composition, while the dashed diagonal 
lines intersecting the solar value are the abundances expected if 
the solar values were scaled with metallicity. 



lation. This enrichment of 12 C results in a very high ratio 
of 12 C/ 13 C. The total carbon, including 12 C and 13 C, shows 
an abundance similar to the observation in Fig. [2] 

Nitrogen is enriched in the great majority of EHe and 
RCrB stars. Heber(1983) and subsequent authors point out 
that the N abundances in general follow the trend of the iron 
abundance. Nitrogen is enriched through CNO cycling in the 
parent stars. In our simulation, it is subsequently reduced 
by the 14 N(a,7) ls O reaction 



Oxygen i s enric hed in most EHe and R CrB stars. |Clay-| 
ton et al. (20071 report that, for four RCrBs, 18 has a 
large enrichment and shows a low ratio of 16 0/ 18 close to 
unity. By comparison with the solar abundance, 18 must 
have increases by a factor > 400. Clayton et al. also indicate 
that the production of 18 O requires temperatures of at least 
10 8 K to allow the 14 N(a,7) 18 reaction. ls O is also pro- 
duced by a-capture on 14 N during the fast-accretion phase 
in our models; the resulting 16 0/ 18 ratio is as a function 
of metallicity, decreasing with Z (Fig. j3|). There is no new 
le O from alpha-capture during the fast-accretion phase. 



Fluorine is 

20061) and ii 



enriched in several EHes (Pandey et al 
most RCrB stars enriched by factors of 



800-8000 re lative to its likely initial abundance (Pandey 
et al. 20081. In our simulation, 19 F comes from reaction 
oT^Nfa, 7) 18 F(p,a) 15 0(a,7) 19 Ne(/3 + ) 19 F. Our simulation 
shows a strong overabundance of F, but not enough to fully 
agree with the observational data. 
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Table 1. Post-merger models with 4.3 ^ logT c ff ^ 4.5. 
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Figure 3. The ratio of 16 0/ 18 as a function of metallicity. The 
computation results are shown as circles with crosses. Diamond, 
triangle and circle symbols show the observations of four 18 O 
enriched RCrB stars , i.e. WX CrA, S Aps, SV Sge and ES Aql 
| |Clayton et al.|2007[ |. 



Neon. A high overabundance of neon has been identified in 
several EHes and RCrBs. In our simulation, 22 Ne is enriched 
by two a-captures on 14 N followed by extensive convective 
mixing. 

Lithium. A few R CrB stars have a notably large overabun- 
dance of lithium, which has so far been difficult to explain. 
There is no significant lithium produced in our simulation. 
|Longland et al.| ( |2012[ ) indicate that lithium can be pro- 
duced by the merging of a helium white dwarf with a carbon- 
oxygen white dwarf if their chemical composition is rich in 
3 He from the previous evolution. This model requires enough 
3 He to be left in the white dwarf after the end of main- 
sequence evolution and a hot enough corona to form during 
the merger. As a test calculation, we put 10 -5 mass fraction 
of 3 He into the accreted material, and obtained a surface 
with a lithium mass fraction of 1.5 x 10~ 5 , or about 10 5 x 
solar. Hence it is clear that if sufficient 3 He is present in the 
accreted helium, He+He white dwarf mergers can also yield 
lithium. Precisely how much depends on many factors, not 
least the quantity of 3 He in the donor and the temperature 
history of the fast-accretion phase. 



3.3 Some Statistics 

The merger frequency of double helium WD systems i n the 

dHan 



Galax y is estimated to lie in 

i 



the range 0.0057 yr 



Han 



(1998) 



1998[ ) (Model 4) to 0.029 yr" 1 ( |Webbink||l984[ ) 
(Fig. 6) shows the final masses of double helium WD mergers 
to lie roughly in a normal distribution with a mean of 0.61 ± 
0.09 Mq. Approximately ~ 2.3% have a merged total mass 
> 0.8 M Q . Thus, the rate for these high-mass mergers should 
lie in the range 1.3 — 6.67 x 10 -4 yr _1 . 

Known low-luminosity helium stars have effective tern- 



Stage 


log L/L 


Timescale (yrs) 


Numbers 


A 


2.7-3 


6 X 10 4 


8-41 


13 


3.6 - 3.7 


8.6 x 10 4 


11 - 57 


C 


4 


2 x 10 4 


3-13 



peratures in the range 4.3 logT e ff ^ 4.5 (Jeffery et al 



19961. Considering the evolution tracks through this tem- 
perature range, we list the luminosity range and evolution 
timescales for three different stages of evolution in Table [l] 
(see also Fig. [T]). Combining the evolution time and esti- 
mates for the merger frequency, we estimate the number of 
low-luminosity helium stars in the Galaxy to be ~ 8 — 41 
in the inwards-burning phase (stage A), ~ 11 — 57 in the 
outwards-burning phase (stage B) and « 3 — 13 cooling to 
the white dwarf phase (stage C). Taking the luminosities into 
account, the relative numbers of massive double- WD merg- 
ers observable as EHe stars in each stage (A:B:C) should be 
~ 1 : 6 : 4, strongly favouring the post-core-burning stages. 

RCrB stars generally have effective temperatures less 
than 10 000K. In the immediate post-merger phase, our 
He+He merger models with T c g < 10 000K have a lifetime of 
~ 6.7 x 10 4 yrs and a luminosity ~ 3 dex above solar. On the 
post-subdwarf giant branch they have an equivalent lifetime 
~ 1.4 x 10 4 yrs and a luminosity ~ 4 dex above solar. Thus, 
we can estimate the Galactic population to be « 9 — 43 new- 
born hydrogen-deficient giants from massive He+He WD 
mergers, and « 2 — 9 more luminous hydrogen-deficient gi- 
ants from post-core-burning He-sdO stars. Assuming no in- 
trinsic extinction, the relative numbers of observable cool 
H-deficient giants in each phase from this channel should be 
» 1 : 3. 

Finally, if we assume a rate of 0.018 — 0.019 yr -1 for the 
merger of CO+He white dwarfs ( [Webbmk|19"84| |Han|199"8] ), 
and a mean evolution timescale which is shorter by a factor 
of 2 due to their higher masses, there would be ~ 14 — 70 
times as many R CrB stars from the CO+He merger channel 
as from the massive He+He channel. 



4 CONCLUSION 

Mergers of helium and carbon-oxygen white dwarfs are cur- 
rently regarded to be the most favoured model for the ori- 
gin of most R CrB and EHe stars. The merger of two he- 
lium white dwarfs has been widely considered responsible 
for the origin of some hot subdwarfs, particulalry those with 
hydrogen-deficient surfaces. In this paper we have demon- 
strated that the most massive helium+helium white dwarf 
mergers can produce hot subdwarfs which subsequently be- 
come cool supergiants and, consequently, provide an alter- 
native means to produce low-mass RCrB and EHe stars. 
However, the number of RCrB stars produced in this way 
may be some 14 to 70 times smaller than from the CO+He 
WD merger channel. Nucleosynthesis of elements during the 
hot (fast-accretion) phase of the merger is roughly consistent 
with the observed abundances of 12 C, 18 0, 19 F and 22 Ne in 
RCrB and EHe stars. Further work is required to estab- 
lish whether the products of massive He+He WD mergers 
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can be clearly distinguished from the products of CO+WD 
mergers. 
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